% SIMPLE 2-D, 1-G, fixed source problem

Y = 10;

it = 1;

%  data format
%    m1
%    g1 -->  D_1  SigR_1  nuSigF_1  xi_1  s_1->s_1 s_2->s_1 ...
%    g2 -->  D_2  SigR_2  nuSigF_2  xi_2  s_1->s_2 s_2->s_2 ...
%    ...
%    m2 ...

dat = [ 1.500   0.030   0.000  1.00     0.000  0.000  
        0.420   0.080   0.135  0.00     0.020  0.000
        2.000   0.040   0.000   0       0.000  0.000
        0.300   0.010   0.000   0       0.040  0.000 ];    

numg = 2;  numm = 2;

xcm  = [   0     2     4      6 ];
xfm  = [     1      1      1    ]*Y;
ycm  = [   0     2     4    6   ];
yfm  = [     1      1     1     ]*Y;
mt   = [  1 1 1 
          1 1 1 
          1 1 1 ];       
src  = [  1 1 1  
          1 1 1  
          1 1 1 ]*0;
     
% xcm  = [   0     2     4      6 ];
% xfm  = [     1      1      1    ]*Y;
% ycm  = [   0     2     4    6   8   10   12 ];
% yfm  = [     1      1     1   1   1    1   ]*Y;
% mt   = [  1 1 1 1 1 1 
%           1 1 1 1 1 1  
%           1 1 1 1 1 1 ];       
% src  = [  1 1 1 1 1 1 
%           1 1 1 1 1 1  
%           1 1 1 1 1 1 ]*0;
      
BCL = 0;
BCR = 0;
BCB = 1;
BCT = 1;           

IBSL = [ 0 ]';
IBSR = [ 0 ]';  
IBSB = [ 0 ]'; 
IBST = [ 0 ]';     